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ABSTRACT 

This paper is the second part of an analysis of the use of a Spacebome Geo- 
dynamics Ranging System for determining crustal strain rates. The present 
analysis focuses on the use of site coordinates rather than intersite base- 
line distances for the strain rate determinations. After discussing the 
analytical techniques which are to be employed, numerical results 
are presented which suggest that the use of site coordinates would result 
in a 20-70$ improvement in the precision of the deduced values of strain- 
ing. Precisions of a few parts in 10 4 would be achievable with simple 
geometries and a decade or two of measurements; precisions of a few 
parts in 10 M would be achievable in a few years. A consideration of pos- 
sible correlations among the derived target site coordinates leads to 
the conclusion that, with the proper choice of coordinate systems, the 
correlations can be made small and non-detrimental to the strain rate 
determinations. 
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DETERMINING CRUSTAL STRAIN RATES 


WITH A SPACEBORNE GEODYNAMICS RANGING SYSTEM, 

2. STATION COORDINATE ANALYSIS 

1. INTRODUCTION 

This paper continues an analysis of the use of the proposed Spaceborne 
Geodynamics Ranging System (SGRS) for determining crustal strain rates. 
This analysis was begun in an earlier paper (Cohen and Cook, 1978) herein 
referred to as Part 1. In l’art 1 we outlined the basic measurement technique 
which would be employed by SGRS The essense of that technique is a very 
precise determination of the locations of a number of targets situated on the 
ground. The survey is conducted by sequentially illuminating the reflective 
targets with a laser pulse from a device located onboard a spacecraft and 
measuring the times of flight for the laser pulses. These time of flights can 
be converted co ranges and with suitable knowledge of the spacecraft orbit, 
the earth's gravity field, and a number of other factors, the targets locations 
on the earth can be deduced. When the grid is resurveyed at a later date, 
changes in the target locations are interpreted in terms of crustal strain 
accumulation. 

In Part 1 we analyzed the precision with which strain rates would be 
deduced from a consideration cf changes in the baseline distances meas- 
ured between the target sites. In the present paper we consider an alter- 
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native method of analysis based on the measurement of the chants in location 
of the sites themselves. This method more fully utilizes the infonnation 
which could be extracted from a SCiKS maasurement. As a consequence the 
deduced formal precisions are somewhat better and additional results, in the 
form of the target grid rotation rates, become available. Conversely, this 
method puts more stringent demands on the accuracy of the SGRS data since 
it depends on the vector resolution of relative Intersite locations rather than the 
single parameter, intersite distance. 

II. ANALYTIC CONSIDERATIONS 

Strain is essentially a geometric concept. It is most easily visualized 
as a distortion of a line between nearby points. This distortion may be in the 
form of a change in the line length or a bending or both. Consider a con- 
tinuous set of points described by the pairs of coordinates (x,y, ) in a two 
dimensional plane. Suppose that we start from an unstrained state and find that 
at some later time a point (x, y) has been displaced to a new location (x-*u,y+v). 
In general the displacements u and v arc functions of the coordinates x and y. 

We can define certain components of an inlinitesmal strain tensor by 


du 

e " = dx 

( la) 

du 

fc * 2 ' = dy 

(id 

to 

II 

( lb) 

dv 

€ 22 ~ * 
dy 

(Id) 

(In Part 1 we used 

tiie symbols e v , e . e 

, ^r e |)t e, 2 . i/ 2 (e 12 + t 21 )). 



It was shown in Part 1 that these components of the strain tensor can be re- 


lated to a change in the distance between nearby points by 


e a = cos J 0 + sin ? 0 + '/j(€, 2 ♦ fj , ) cos 0 sin 0 (2) 

r o 

whore r is the distance between the points In the strained state, r a , is the 
unstrained distance, and 0 is the angle the interlocution line makes with the 
positive x axis. 

In general the unstrained state of the surveyed region is not known; 
rather the region s surveyed at different times to determine the rates of 
strain aecumulatio.i. Differentiating equation 2 we find 

r = r %1 j cos J 0 + e , 2 sin* 0 + 'jir, , + ) sin 0 cos 0 1 (3) 

Equation 3 suggests one method for determining strain rates. By determining 
changes in the lengths o* three or more appropriately chosen lines in a region 
of homogeneous straining, the components of the two dimensional strain rate 
tensor can be deduced. We showed in Part 1 that with relatively simple site 
configurations, modest baselines of 25-70 km., and typical baseline pre- 
cisions of a few centimeters, the strain rates could be determined to a pre- 
cision of several parts in 10 v in a measurement program lasting several 
years. 

There are alternative procedures for deriving strain rates depending 
on the survey data available. In the present analysis we will assume that 
each of the target locations are known from the survey data in some conven- 
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tent reference frame: some discussion of the selection of the reference 
frame Is presented later In this paper. Consider the x coordinate of two 
points ( x, , y,) and x 2 , v,). Assuming that the coordinate change linearly 
with time 

x, <t + At> - x,(t) * u(x r y,)At (4) 


Furthermore to lowest order 


aii aii 

An = u(x,.y 2 ) - ut\, ,y, » “ — <x 2 - + — <y> - Yj > 


ay 


( 5 ) 


i. e. , 


Similarly 


Aii = e,, Ax + f,, iy 


Av = e, 2 Ax + c, 2 Ay 


<<>> 


(7) 


The preceding equations form the basis for the strain-rate analysis 
using station coordinates. Suppose as a result of several resurveys one has 
a collection of coordinates (x^), y^t-)) for several target locations at 
several times tj then it is a straight forward matter to calculate u and v. For 
example, using least square techniques if 


xttji = a + lit. 


(S) 
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then for n measurements 
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Given the rates u and v we can fo**m differences Au,. Av ( for the rates between 
pairs of sites. Then if there are p pairs of target sites 


n 




, 11,1 


£au Ax. Av Ax 
, i i , i i 

lAu, Ay, b», Ay, 


(10) 


This equation should lx; discussed in some detail The left hand side contains 
four components of " tencrulized strain rate tensor. We can write this tensor 
as a sum o * a symmetric and antisymmetric tensor. 
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The symmetric tensor is the pure strain rate tensor discussed in Part 1, 
while the antisymmetric tensor represents a rotation rate. 

Next we note that the inverse matrix on the left side of equation 10 
gives the geometric contribution to the variance-covariance matrix for the 
strain rate determination When multiplied by the variance in the coordinate 
rate differences, the diagonal terms give the variances in the strain rates and 
the off-diagonal terms give the covariances between terms. Let 
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then the inverse matrix is (A *A) 1 and 
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where the r, are interstte distances and the 8, and m, are direction cosines 
with respect to the x and y axes respectively. Let o*. and o‘ t be the variances 
in au and Av and o‘ x and o 2 y be the corresponding variances in Ax and Ay. 

From equation 9 N > 

n , 


°lu 


nit? - (it > 2 

l 1 ( ’ 


AX 


( 14) 


The variances in the strain rates can be calculated from equations 13 and 14, 
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Furthermore If we assume and t ,, are uncorrelated the variance in the 
shear strain is 1 .•(<>: ♦ oj ) which is also the variance in the rotation rate. 

' i: ' :i 

In the next section we will present results showing the expected standard 
deviations in the strain rate dele* mlnaHons. These values are simply the 
square roots of the preceding variances. 

III. RESULTS 

The precision with which the strain rates can be deduced depends on the 
geometry of the target grid, the number of sites, the resurvey period, At, 
the period between the first and last measurements, T, and, of course, the 
precision In the deduced site locations. In Figure 1 we show five simple 
target grids for which wo have calculated the expected precision using the 
present analv^ls and which we also used in Part 1 for that analysis. The 
results of our calculations are shown in Figures 2-6. These results have 
been normalized to a one centimeter standard deviation in the difference in site 
coordinates, a normalization consistent with the one centimeter standard deviation 
in baseline distances used in Part 1. (If we assume that all coordinates are 
determined to the same precision and are uncorrelated then, for example, 
o„- a, - y/f o % .) For the simple triangle of Figure la the precision in the de- 
duced strain rates improves from about 2 parts in 10 7 per year after tw’o years to 
1-3 parts in 10 K per year after 15 years. These numbers should be compared 
to the typieal rates of straining on, say, the San Andreas Fault of a few parts 
In 10 7 per year. Considerable improvements in the precision can be obtained 


In going to ;;ri«N Involving more than three targets For the 9 target grid of 
Figure le, the precision varies from about 4 parts tn 10* per vear after 
two years to 2-3 parts in 10* after fifteen years 

It is interesting to compire the results obtained hen* using site eoordtnntes 
with thole obtained from a baseline distance analysis, -jell a comparison is 
shojvn In Table 1 The principal strain rates, e,, and c ,,, derived rom the 
coordinates are typically 20 percent more precise than those obtained using 
baselines. The shear strain rate precisions Improve by a factor 1.4- 1.7 in 
going from the baseline to coordinate analysis. Furthermore the relative 
precision in the deduced values ot the components of the strain rate tensor 
mat vary with the method of analysis. For Figure Id the precision in 
is twice that of t‘,, in the site coordinate analysis butv^» better in the baseline 
distance analysis 

We have assumed in the analysis presented here and in Part 1 that the 
correlations in the deduction of various components of several site locations 
are negligible. Wo should consider this paint in a bit mo-e detail. If positive 
correlations exist among common components, sav the x component, of 
different station locations, then Uv* measurement and analysis errors give 
rise to a common po sition bias which is largely removed when the differences 
between target locations arc determined in cith ; o*. the analysis methods we 
have considered. On the other hand if the correlations are large and negative, 
the precisions can be substantially degraded. We examine the question of 
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correlation* bv considering a varlance-co variance matrix P for the site 
location*. Let 

Ax 

h 

Ax 

S ■ Ay 

• 

• 

Ax 
Ay 

be the vector of uncertainties in the site locations. The matrix P la defined 
by 

P * <ss r > 

( 17 ) 

Associated with the variance-covariance matrix is a correlation matrix, 
with definition 



= P (P P 

I) II D ' 


The diagonal elements of the correlation ma< rlx have a value of unity while 
the values of the off-diagonal terms vary from +1 (fully correlated) to -1 
(fuliy correlated negatively) with the value 0 corresponding to no correlation. 


If the quantities a and & have a correlation p o(J then 

o- = o* + or ± Zpoa.. 

atfi a J a iS 


In order to develop a quantitative appreciation for the correlations we 


consider a network of nine stations which are part of a larger grid being 
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considered for part of California. The coordinates for the targets are shown 
in Table 2; the arrangement of sites is similar to that of Figure le although 
the sides of the grid are not along East-West and North-South lines. The 
expected precisions in site locations and the correlations among coordinates 
have been considered in two separate coordinate systems. The first system 
is an earth center fixed coordinate system (ECF) in which the z axis is aligned 
with the Earth’s rotation axis, the x axis pointa toward the prime meridian through 
Greenwich and the y axis completes a right handed orthogonal system. The second 
system is a local tangent plane coordinate system (LTP) with one of the target 
sites held fixed at the origin. In this system the x axis points to the East and the 
y axis to the North. The z axis is vertical. The mean values of the 
correlations determined in each of these two systems are shown in Table 3. 

(Not shown are the r.m.s precisions in the position determinations which are 
typically 11 cm. for each component in the ECF system and about 2 cm. for 
each component in the LTF system.) On comparing these results we conclude 
that the balk of the errors in the site locations in the ECF system are bias 
errors which are partially removed by transforming to the LTP system. 

In the LTP system mean values of the correlations are less than 0.3 except in 
the vertical direction. It may be significant to note that the correlations among 
the x and y components have nearly zero means while the correlations between 
common components of different stations still have a positive mean. This suggest 
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the presence of a residual bias which might be removed by transforming to a 
new coordinate system with a resultant further Improvement In the precision. 
For the purposes of this paper, however, the correlations are already suffi- 
ciently small to suggest that our deductions concerning strain rate precisions 
should be correct to with a factor of two or so. 

IV CONCLUSIONS 

We have considered the potential use of measurements made with a SGRS 
system for determining crustal strain rates. In Part 1 a method of analysis 
was outlined which focused on the use intersitc distance measurements while 
In Part 2 the method relied on a decomposition of the Intersite vector Into 
relative coordinates. Numerical calculations using small target grids with 
25-70 km site separations suggest that straining rates can be determined to a 
few parts in 10 8 per year, after several years. Precisions of several parts 
in 10 v may be attainable with a decade or so of measurements. The use of 
station coordinates as opposed to baseline distances result in a modest 
improvement In the precision of the derived strain rates as well as permitting 
more detailed decomposition of the strain rate tensor. 
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TABLE 1 

— 

Relative Standard Deviations in Strain Rate Determinations Based oa Relative 

Site Coordinates, 

o c , and Intersite Baseline Distances, 

V 



Component 


Figure no. 

Jn 

«22 

*<«I2 + «2I> 

_ 

la 

G b / o c =1. 2 

0 b /0 c «1. 2 

o b /o c «L 7 

lb 

1.2 

1.2 

1.7 

lc 

1.2 

1.2 

1.4 

Id 

1.1 

1. 3 

1.4 

le 

1.2 

1.2 

1.5 
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-TXBEET 

Station number 

Longitude 

Latitude 

Elevation (m) 

151 

119*57’ 

36*01* 

61 

152 

119*40 * 

35*50' 

304 

153 

119*36* 

35*40' 

910 

161 

120*10' 

35*52’ 

305 

162 

119*59' 

35*42' 

215 

163 

119*48' 

35*31' 

152 

171 

120*23' 

35*43' 

457 

172 

120*11' 

35*32' 

609 

173 

120*00' 

35*22' 

610 



TABLE 3 


Correlation coefficients associated with off diagonal terms of coordinate 
covariance matrix calculated for a nine station subnetwork of proposed 
California grids. Station 162 is the fixed station for calculations in the 
I TP system. 

ECF Coordinate System LTP Coordinate System 
Compone nt <p) <{p-<p> i )> Vi <p> <(p-<p>) 2 >* 


V-I ] 

.97 

CM 

O 

• 

a 7 

.37 

y , y , 

.98 

.01 

.29 

.32 

Vj 

.99 

\ i *j 

.01 

.48 

.01 

Vj 

.06 

.08 

.02 

.35 

X i*j 

.08 

.09 

.09 

.02 

y. z i • 

.36 

.04 

- k 09 

.02 

v. 

.06 

.07 

-.00 

.24 

Vi 

.09 

.09 

.19 

.08 

y. z . 

.36 

.06 

-.17 

.04 
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FIGURE CAPTIONS 


Figure 1. 
Figure 2. 

Figure 3. 

Figure 4. 

Figure 5. 

Figure 6. 


Representative SC1RS target site grids. 

Calculated precision in the strain rates c, , , e 2| , e |2 , «* 22 for the 
grid shown in Figure la. 

Calculated precision in the strain rales f,, , t,, ,e, 2 / 2, 1 22 / 2 for 
the grid shown in Figure lb. 

Calculated precision in the strain rates i , , , e* 2 , , c, 2 , « 22 for 
the grid shown in Figure lc. 

Calculated precision in the strain rates r, , , e-,, , ^, 2 /2,e, 2 /2 for 
the grid shown in Figure Id. 

Calculated precision in the strain rates , , t' 2| , e, 2 ,e 22 f° r 


the grid shown in Figure le. 



e 

Figure 1 Representative SGRS target site grids 
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Figure 3. Calculated precision in the strain rates 
e,,» e 21 , e, 2 /2, e 22 /2 for the grid shown in Figure lb 
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T lyMitl 

Figure 5. Calculated precision in the strain rates 
C| | • c 2 , , for the grid shown in Figure Id 
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Figure 6. Calculated precision in the strain rates 
e n , e 21 , e. 2 » e 22 for the grid shown in Figure le 



